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ABSTRACT 

Intermediate resolution spectra at the Ca II triplet have been obtained for 55 candi- 
date red giants in the field of the globular cluster M22 with the VLT/FORS2 instrument. 
Spectra were also obtained for a number of red giants in standard globular clusters to 
provide a calibration of the observed line strengths with overall abundance [Fe/H]. For 
the 41 M22 member stars that lie within the V — Vhb bounds of the calibration, we 
find an abundance distribution that is substantially broader than that expected from 
the observed errors alone. We argue that this broad distribution cannot be the result 
of differential reddening. Instead we conclude that, as has long been suspected, M22 
is similar to oj Cen in having an intrinsic dispersion in heavy element abundance. The 
observed M22 abundance distribution rises sharply to a peak at [Fe/H] pa -1.9 with a 
broad tail to higher abundances: the highest abundance star in our sample has [Fe/H] 
pa -1.45 dex. If the unusual properties of lo Cen have their origin in a scenario in which 
the cluster is the remnant nucleus of a disrupted dwarf galaxy, then such a scenario 
likely also applies to M22. 

Subject headings: globular clusters: general — globular clusters: individual (M22) - 
stars: abundances 



1. Introduction 

Beginning from observations made more that 30 years ago, we now know that the majority, if 
not all, globular clusters show star-to-star variatio ns in the abundances of the light elements C, N, 



O, Na, Al and Mg (see, for example, the reviews of lKraftlll994l ; iGratton et al .1120041 . and references 
therein). The variations are in the sense that relative to 'normal' stars in the cluster, carbon and 
oxygen are depleted while nitrogen is enhanced, and sodium and aluminum are enhanced while 
magnesium is depleted, in the stars showing the 'anomalous' abundances. Together these effects 
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are known as the O-Na anti-correlation. Despite much work, the origin of the abundance anomalies 
is not well understood though the fact that they are seen on the main sequence in at least some 
globular clus ters points to a proce ss that is intrinsic to the formation of the cluster (e.g., see the 
discussion in lMarcolini et al.ll2009l ). 



In general, however, globular clusters are chemically homogeneous when it comes to the abun- 
dances of the iron-peak elements, with the limits on any possible internal range in [Fe/H] quite 
stringent. For example, in their detaile d analysis of h i gh di spersion spectra of 36, 23 and 28 red 
giants in the clusters M5, M3 and M13, iKraft fe Ivansl (|2003l ) list observed standard deviations for 
the cluster [Fe/H] values in the range 0.03-0.08 dex, consistent with that expected from the errors 
alone. Consequently, any intr insic abundance spreads must be considerably smaller. Similarly, 
Da Costa fe Armandrofa ()1990l ) used their red giant branch photometry to set 3a upper limits of 
0.04-0.09 dex for any intrinsic heavy element abundance range in six southern clusters. 

The well-established exception to the lack of [Fe /H] variations in individual clusters is the stellar 
system u Centauri, the most luminou s of the Galaxy's globular clusters. While the stars in this 
cluster show the O-Na a nti-correlation (INorris fe Da Costalll995aT). they also possess a wide range in 
overall abundance (e.g.. iFreeman fe Rodgerslll975l ; INorris fe Da Costalll9 95b; Pa ncino et ajj [2002j) 



Romano et al 



2007, and 



with complex distributions of element-to-iron abundance ratios (see 
the references therein). This indicates that uo Cen is a system that has undergone significant 
chemical evolution. Indeed the nucleosynthetic contributions of Type II and Type la supernovae 
are recognizable in the variations of abundance ratios with [Fe/H], as are the c ontributions of 



AGB stars (e.g.. INorris fe Da Costal Il995bl ; ISmith et alJ 120001 ; IPancino et alJ 12002 ). All in all the 



observations point to a complex chemical history for uj Ce n, with star formati on likely to have 
occurred in the cluster over a period of perhaps 2 Gyr (e.g., IStanford et al.l 12006). 



The differences between to Cen and other globular clusters has led to the suggestion that u Cen 
may have formed in a different way — that it is the nu clear remnant o f a dwarf galaxy which ha s 



been disrupted by the tidal field of the Milky Way (e.g., lFreemanlll993l ). iBekki fe Freeman! ([2003), 
for example, have shown that such a process is dynamically plausible. The different environment 
may then ha ye facilitated additional chemical processe s that do not occur in 'regular' globular 
clusters (e.g.. lBekki fe Norrisll2006l ; iRomano et al.ll2007l ). 



In this context the location of the globular cluster M54 at the center of the Sagittarius dwarf 
spheroidal galaxy is particularly relevant. The Sgr dwarf is currently being tidally disrupted by 
the Galaxy and in a few Gyr or less, M54 will be seen as a halo globular cluster rather than as the 
central star cluster of a dwarf galaxy (e.g., Da Costa fe Armandrofjll995 ). It is therefore potentially 
a current day example of what may have been the situation for uj Cen in the distant past. From 
an abundance point-of-view the situation is complicated because the cluster is superposed o n the 
general Sgr field population, as well as on that of the Sgr nucleus. iBellazzini et al.l (|2008l ) have 
shown that the Sgr nuclear population is metal-rich ([Fe/ H] ~ -0.4) and kinem a tically distinct from 
that of M54, which has [Fe/H] ~ —1.5 dex. Nevertheless ISaraiedini fe Lavdenl (|l995l ) claimed from 
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their analysis of the width of the M54 red giant branch in a (I, V — I) co lor-magnitude diagram that 
cr([Fe/H])j n t « 0.16 dex in M54. More recently iBellazzini et al.1 (j2008l ) have obtained Ca II triplet 
spectroscopy for over 700 red giants in the central M54/Sgr field. They associate ~425 stars with 
M54 and find <r([Fe/H])j n t ~ 0.14, in good agreement with previous estimates, and an observed 
range in [Fe/H] from -1.8 to —1.1 dex. Thus the case for M54 being the second star cluster after 
<jj Cen to possess an internal range in [Fe/H] is strong, but confirmation from an extensive high 
dispersion spectroscopic study remains to be done. Such a study would also allow measurement of 
[element/Fe] ratios and their variation (if any) with [Fe/H], from which, as for uj Cen, constraints 
could be placed on the enrichment processes. 

The other globular cluster that is often mentioned in the context of uj Cen-like abundance 
variations, albeit on a smaller scale, is M22 (NGC 6656). Unlike uj Cen and M54 which are 
both very luminous clusters with My ~ -10.3 and -10.0, respe ctively, M22 is a cluster of bright 
but not outstanding luminosity, with My ~ -8.5 (jHarrid I1996I ). It lies at low galactic latitude 
approximately 3 kpc from the Sun and 5 kp c from the Galactic Center ([Harris! 11996). Based 
on DDO photometry for 10 M22 red giants, Hesser. Hartwick McClure ( 19771 ) were the first 
to suggest an analogy between M22 and uj Cen by noting that the observed range of ultraviolet 
excesses and CN-stren gths in their M22 sample was similar to that seen in uj Cen stars (see also 
Hess er fe Harrislll979l ). This was followed by the much more extensive study of lNorris &i Freeman 
(|1983l ). who obtained low resolution spectra for ~400 M22 red giants. They reported the existence 
of a correlation between the strength of cyanogen features and the strength of the Ca II H and K 
lines in their M22 spectra: such a phenomenon was also seen in similar spectra of u Cen red giants 
but wa s much less marked for r ed giants in the 'normal' cluster NGC 6752. Based on synthetic 
spectra iNorr is Freeman! (|1983l ) then went on to conclude that the observed range in calcium line 
strengths in M22 corresponded to an abundance range A[Ca/H] ~ 0.3, and that the cluster shared 
the anomalous abundance patterns of ui Cen, though to a smaller degree. 

In the ensuing decades the debate as to whether M22 does or does not share many of the 
characteristics of uj Cen, particularly as regards the presence of a range in heavy element abundance, 
has swung back and forth with little consensus and a number of divergent results. The situation is 
complicated by the clea r presence of differential reddening across the field of the cluster: AE(B — V) 
~ 0.06 - 0.0 8 mag (e.g., Anthony-Twarog et alll9"95 : Monaco et al. 2004 ). A brief sur vey of existing 



result s (see iMonaco et al.1 |2004j, for additional references) includes the following. iLehnert et al 



(|199ll ) used spectra at Ca II triplet region of 10 red giants to conclude that M22 was similar to uj Cen 
in displaying Ca, Na and F e abundance variat i ons. T hey estimated A[Ca/ H] ~ A [Fe/H] ~ 0.4 dex, 
a resu lt similar to that of lNorris fc Freeman! (|1983l ). On the other hand, lAnthonv-Twarog et al 



(|1995l ) concluded from Stromgren+Ca photometry of ~300 r ed giants and horizon tal branch stars 
that there was no evidence for a range in [Fe/H] in the cluster. IMonaco et al.l (|2004l ) reached similar 
conclusions from their extensive wide-field photometr y of M22: the maxim um metall icity spread 
permi tted by their data is A [Fe/H] ~ 0.1 - 0.2 dex ([Monaco et al.ll2004l ). Further, llvans et al 



(J2004J), using high dispersion spectra for 26 M22 red giants, were able to find a consistent set of 
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spectroscopic and chemical constraints which gave acceptable stellar parameters and no requirement 
for any variations in [Fen/H]. 



Most recently, iMarino et al.l (|2009l ) analyzed high dispersion UVES spectra for 17 red giants, 
and lower resolution GIRAFFE spectra for 14 stars, in the cluster. In addition to the O-Na anti- 
correlation that is also seen in many clusters, they identified the presence of two groups of stars 
whose mean abundances for the s-process elements Y, Zr and Ba differ by ~0.6 dex. The s-process 
rich group also ap pears to have highe r iron and calcium abundances by 0.14 ± 0. 03 and 0.25 ± 0.04 
dex, respectively (jMarino et al.l l2009). In this context it is also worth noting that iPiottol ((2009) has 
presented a HST-based color-magnitude diagram for M22 in which it is evident that there are two 
distinct cluster sub-giant branches. M22 thus joins other clusters such as u Cen, NGC 2808 and 
NGC 1851 where such diagrams provide evide nce for the p resence of two or more internal cluster 



populations that have different properties (see |Piottoll2009l . and references therein 



As part of a larger program to study Galactic globular clusters with uncertain abundance deter- 
minations, we have investigated the question of an abundance spread in M22 anew, by obtaining in- 
termediate resolution spectra at the Ca II triplet of a substantial number of M22 red giant stars. Ca 
II triplet line strength measu rement from such spectra is a well established technique for determin- 



ing o verall abundances (e.g.. lArmandroff &; Da Costa 



1991 



Rutledge et al.l ll997b; B attaglia et al 



20081 ) . In §2 the observations and reductions are described, while in §3 a calibration of the measured 



line strength indices with overall abundance is generated. The results for M22 are also presented 
in §3 and are discussed in §4. In brief, we find strong evidence that there is an overall abund ance 
spread in our M22 s ample of size similar to that originall y found by iNorris Freeman! (119831 ) and 
Lehnert et al. ( 1991 ). and consistent with the results of Marino et al. ( 20091 ). Moreover, the M22 



abundance distribution bears some similarity to that for u Cen, although on a reduced scale. 



Observations and Reductions 



2.1. Observations 



Short and 'long' exposures of each cluster in the program were first obtained in service time 
with the FORS2 instrument in imaging mode at the Cassegrain focus of the VLT-UT1 telescope, 
primarily to define the targets for spectroscopic follow-up. The 6.8' x 6.8' field-of-view was centered 
on the cluster for the sparser or more distant systems, but was offset from the centers for the 
nearer and richer clusters, including the 'standard' clusters. These latter clusters w ere taken from 
the c ompilation of lPritzl et al.1 (|2005l ) and possess well established [Fe/H] values (e.g. lKraft k, Ivans 
2003) that cover the full range of metallicities exhibited by Galactic globular clusters. Observation 
of such clusters in addition to the program clusters allows the derivation of a line strength - 
abundance calibration from the subsequent spectroscopic observations. The images were obtained 
in the V and I bands. 
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Point-Spread-Function photometry was carri ed out on both the long and short exposure pairs 



using Stetson's DAOPHOT / ALLSTAR package (|Stetsonlll987l . Il994h . The resulting instrumental 



magnitudes were then provisionally calibrated by using color terms, which are small, and zero 
points provided by ESO as part of their routine quality control. Color-magnitude (c-m) diagrams 
were then generated to allow the selection of likely cluster members as targets for the spectroscopic 
observations. They also permit characterization of the principal c-m diagram features, such as 
horizontal-branch morphology. We note though that for the clus ters considered here , superior 



photometric data are available from other sources, for example from lMonaco et al.1 (120041 ) for M22. 
The calibration approach adopted means that the photometry is likely to be consistent with the 
standard V, I system, but that the zero points for each cluster data set are uncertain at the 0.05- 
0.10 mag level. However, since our analysis uses only differential magnitudes, such as the magnitude 
difference from the horizontal branch, the zero point uncertainty is not a concern. We note that 
for each cluster the data sets from the two detectors in the camera were treated separately: stars 
falling on CCD1 are labeled as '1_' followed by a running number and those falling on CCD2 as '2J 
(cf. Tables EE}. 

The spectroscopic observations were then carried out in visitor mode with the FORS2 instru- 
ment in MOS-mode during a two night run at the end of May in 2006. Conditions on the first 
night were photometric while on the second night there was intermittent thick cloud. On both 
nights the seeing was quite variable, from below 1" to worse than 2". The instrument was used 
with the 1028z+29 grism and the OG590+32 order-blocking filter. This gives a maximum spectral 
coverage of ~7700A-9500A at a scale of 0.85A per (binned) pixel. The allocation of targets to slits 
was done in such a way that the wavelength interval ~8200A-8900A was always on the detector. 
The number of stars observed in each cluster per configuration varied from ~10 to the maximum 
of 19, depending on the density of likely members within the field-of-view. The magnitude range 
was set as ~3 mag (or somewhat less) to insure sufficient coverage of the red giant branch, and 
the exposure times were chosen to ensure the brightest stars were close to, but not saturated. Two 
exposures were obtained to allow removal of cosmic-rays. For the majority of clusters only single 
configurations were used, but for M22 three separate configurations were observed to increase the 
sample size. The M22 field was centered to the east of the cluster center (see Fig. [5|). 



2.2. Reductions 



All spectra were extracted using the FORS2 pipeline version 1.2 (jlzzo &; Larsenl l2008). and in 
particular the FORS_CALlB and FORS^SCIENCE pipeline recipes were used to wavelength calibrate 
and extract the spectra. The first recipe takes as input the master bias frame, the screen flat-field 
frames and the arc frame, and computes the wavelength calibration and distortion map for each 
slit, after finding the slit positions. A catalog of arc lines and a grism parameter table, part of the 



2 see www.eso.org/observing/dfo /quality /FORS2/qc / photcoeff /photcoeffs_fors2.html 
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package, are also given as input. The mean residual in the wavelength calibration was 0.24 pixels, 
the mean spectral resolution was R ~ 2440, and the mean fwhm of the arc lines was 3.51±0.07 A. 

Once the calibration tables were created, they were used with the forsjscience recipe to 
reduce the target spectra. This recipe cannot process multiple exposures simultaneously, so each 
frame was reduced independently, and then the two exposures were average-combined. The software 
corrects for bias and flatfield, and computes a local sky background for each slit, to be subtracted 
from the object spectra. Cosmic-ray rejec tion was not applied. The extraction radius was 6 pix, 



and Home optimal extraction was applied (|Horndll986l ). The spectra were normalized by exposure 



time, and the wavelength solution was aligned to a reference set of >20 sky lines by applying 
an offset. The median offset was 0.9 pix. The final product of the pipeline is then a fits cube 
containing the fully reduced spectra. Each object is identified by its position and extraction window 
in the rectified spectrum, and by the row number in the fits cube holding the reduced spectra. 
To associate each spectrum to the corresponding target star, the spatial transformation between 
the 'rectified CCD' coordinates of the slitlets and the original target coordinates were computed. 
In this way it was easy to verify that the objects detected in each slitlet corresponded to the 
expected target. Finally a table was created, where for each star the J2000 coordinates, the V and 
/ magnitudes, and the name of the associated spectrum are given. 

The S/N ratios for the final spectra varied from ~110 for the brightest stars to ~25 for the 
faintest in a typical exposure. In Fig.[T]we show examples of the final spectra for a star in NGC 6397, 
two stars in M22, and one star in M10. All four stars have V— Yhb ~ ~0.5 mag. 



3. Analysis 



The first step in the analysis of the reduced spectra for each cluster was measurement of the 
individual radial velocities. This was carried out with the rvidlines task in IRAF's RV package 
using all three of the Ca II triplet lines. After heliocentric correcti on the i ndivid ual velocities were 
compared with the cluster radial velocity from the compilation of lHarrisI (|l996l ). Obvious outliers 
were discarded as definite non-members and the remaining velocities averaged. For the 5 standard 
clusters (M15, NGC 6397, M10, M4 an d M71) , the m ean difference between the cluster velocities 
determined here and that tabulated by lHarrisI (1996) was ± 3 km s -1 , where the uncertainty is 
the standard error of the mean. 



3.1. Abundance Calibration 



For the radial velocity member stars in the standard clusters the (pseudo) equivalent widths 
of the A8542A and A8662A lines of the Ca II triplet were measured by fitting gau ssian line profiles 
using the feature and continuum bands listed in lArmandroff &: Da Costal (| 19911 ) . Along with the 
value, the fit also returns an uncertainty in the measured equivalent width. The sum of the 
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Fig. 1. — Example spectra from the clusters NGC6397, M22 and M10. All 4 stars have approx- 
imately the same V-Vhb value. The spectra have been normalized at A w 8600A but have not 
otherwise been altered from those used in the analysis. 
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two measurements, W8542 + ^8662, was then plotted against the magnitude difference from the 
horizontal branch, V — Vhb, with the measurement uncertainty in Wg542 + ^8662 taken as the 
quadratic combination of the uncertainties for the two lines. For all clusters the V values come 
directly from the photometry of the FORS2 imaging observations with the Vhb values determined 
from the color-magnitude diagrams generated with those photometry sets. As noted above, the 
CCD1 and CCD2 sets were treated separately though in practice the two Vhb values generally 
agreed well. Three M71 stars, and one NGC 6397 star, lay off the (Ws542 + Ws662> V — Vhb) 
relations defined by the other cluster stars; we assumed these stars are non-members. Given that 
both M71 and NGC 6397 have relatively low radial velocities, the fact that velocity is an imperfect 
membership discriminator is not surprising. 

To increase the number of standard clusters available to establish the abundance calibration, 



the standard cluster observations from lGullieuszik et al.l (|2009l ) were also analysed. These spectra 
were obtained with the identical instrument setup and were reduced using the same methods as 
employed here. The clusters observed were NGC 4590, NGC 4372, NGC 6397 (no stars in common 
with our data set), NGC 6752, M5 and NGC 6171. The line strengths on these spectra were 
measured in an identical fashion to those for the standard clusters observed here. The V — Vhb 



values were taken from lGullieuszik et al.l (|2009i ). 



For each of the 10 standard clusters, the ( Ws542 + ^8662 > V — Vhb) data were fit with a straight 
line via least squares. The dependence of the resulting slopes on abundance was investigated and 
found to be negligible. A weighted average of the slopes was then formed and the resulting single 
value of -0.51 ± 0.01 A/ mag refitted to the data for each cluster. The results of this process are 
shown in Fig. [2j Two points are worth noting. First, as is visible in Fig. [2] in the data for NGC 
6397 and to a lesser extent for M10, the slope of the (Ws542 + W'8662, V — Vhb) relation appears 
to notably flatten for V — Vhb values greater than approximately +0.2 mag. Consequently, in the 
analysis here and subsequently, only stars with V — Vhb less than +0.2 are considered. Second, 
the slope derived here is larger than the values near -0.64A /mag found in other studies (e.g . , 



Armandroff & Da Costa 1991 



Rutledge et al.lll997al : iBattadia et al.ll2008l : iGullieuszik et al.ll200fll ) 



Our explanation for this difference is the following. Because o f the s omewhat higher resolution of 
these spectra compared to those of lArmandroff fc Da Costal (119911 ) for example, the gaussian- 
only fitting technique employed here underestimates the contribution to the pseudo-equivalent 
widths from the line profile wings, with this underestimate being larger for stronger lines (see 
Battaglia et al.l 120081 . for similar comments) . We have in yestigated this th r ough a comparison 
of the Ws542 + W8662 values measured here with those of IGullieuszik et al.l (|2009l ) , who used a 
'gaussian + lorentzian' measurement technique that includes the contribution from the line wings. 
The comparison shows an excellent correlation over the full range of measured value s for th e 41 stars 
in 6 clusters in common. However, the Wg542 + Ws662 values of IGullieuszik et al.l (120091) are 17% 



strong er, and when coupled with the 0.51 A/mag slope found here, predicts that 



Gullieuszik et al 



(200 9j) should see a slop e of 0.60 A/mag for their relation between Ws542 + ^8662 and V — Vhb- 
Gullieuszik et al.1 ( 20091 ) actually determine a slope of 0.63 ± 0.02 A/mag. Thus our explanation 
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for the smaller slope found here is likely the correct one. However, since there is no dependence of 
the fitted slope on abundance, and since the range of V — Vhb values covered is similar for most of 
the clusters, our use of the slope determined here should not affect any of the results. In particular, 
we note that if we consider only stars in the standard clusters with V — Vhb < -0.5, rather than 
< +0.2, then the weighted mean slope found is -0.54 ± 0.02 A/ mag. This value is not significantly 
different from that adopted and indicates that our choice of a limit at V — Vhb =0.2 for stars to 
include in the abundance calibration process does not unduly influence the results. 

The variation of Ws542 + 1^8662 with V — Vhb can then be removed by calculating the value 
of the reduced equivalent width W , where W is the mean value of W8542 + W8662+0-51(V — Vhb) 
for the stars in each cluster, subject to the condition that V — Vhb < 0.2 where necessary. This 
single parameter W' then needs to be calibrated with some other measure of cluster abundance. 
In most situations the calibration is made to overall abundance, designated by [Fe/H], with the 
[Fe/H] yalues chosen f rom o ther compilations. Armandroff Da Costal ( 1991 ) ado pted the [Fe/H 



scale o f IZinn Sz Westl (119841 ) while others, for example IRutledge et al 



(1997 bl) and 



20081) . have provided calibrations to the [Fe/H] scale of Carretta Gratton (|1997 ). 



Battaglia et al 



Kraft Sz Ivan; 



2003) have derived a new [Fe/H] scale for globular cluster abundances by consistently analyzing 
high-dispersion spec tra for a number of red giants in 16 key clusters with -2.4 < [Fe/H] < -0.7 



dex. In their study iKraft &: Ivansl (|2003l ) consider the effects of different model atmospheres and 
different color-temperature relations. Here we adopt their [Fe /H] values bas e d on the MARCS 
model atmospheres as our fundamental calibration points. As IKraft 8z Ivans! (12 003) note, these 
abundances are generally about 0.2 dex lower than those of ICarretta &: Grattonl (119971 ). 



Fig. [3]shows the relation between W' and [Fe/H] ma^cs from lKraft fc Ivansl (|2003l ) for our 10 
standard clusters. Over the range of [Fe/H] MARCS shown, the relation is linear with no indication of 
any change in slope towards the higher metallicities. We defer to a subsequent paper the question of 
the calibration for clusters that exceed the abundance of M71, the most metal-rich cluster included 
here. The error bars for the W values come from the dispersion of the H / 8542 + ^8662 values about 
the fitted slope (cf. Fig. [2]) while the erro rs in the [Fe/H]MARCS values are taken from the discussion 
in the Appendix of lKraft Ivansl (|2003l ). A linear regression applied to the points in Fig. [3] yields 
the equation: 

[Fe/H] marcs = (0.549 ± 0.031) X W' - (3.369 ± 0.102) (1) 

The rms dispersion about th e fitted line is 0.092 d ex which is consistent with the average uncertainty 
of the [Fe/H] marcs values (jKraft &: Ivansl 120031 ). and we adopt this relation as our calibration to 
overall abundance. To within the uncertainties our calibration relation is identical to that, namely: 



[Fe/H] marcs = (0.531 ± 0.025) xW- (3.279 ± 0. 



(2) 



given b y lKraft &: Ivansl (120031 ) for the relation between [Fe/H] a/ arcs and the W values of lRutledge et al 
(|l997ah . 



The data for the standard cluster stars observed here are given in Table [TJ In successive columns 
we list the cluster and star ID, the RA and dec, the heliocentric velocity, the value of V — Vhb arid 
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Standard Cluster Data 



«1r 




M4 

NGC6171 
M5 



NGC6752 
M10 



NGC6397 
NGC4372 



M15 

NGC4590 



Fig. 2. — Plot of Ca II line strength (H / 8542+M / 8662) against magnitude difference from the horizontal 
branch (V — Vrb) for the 10 standard clusters. In order of increasing Ws542 + Ws662 values at 
V — Vrb = -1.5, the solid lines are for clusters NGC 4590 (individual stars plotted as open 6-pt 
star symbols), M15 (open triangles), NGC 4372 (filled triangles), NGC 6397 (open 5-pt stars), M10 
(open circles), NGC 6752 (filled 5-pt stars), M5 (plus symbols), NGC 6171 (open squares), M4 
(x-symbols), and M71 (open diamonds). The data for each cluster has been fit with a line of slope 
-0.5lA/mag for V — Vhb < 0.2. Vertical bars on each point show the measurement uncertainty in 
the line strengths. 
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Fig. 3. — Plot of the reduc ed equivalent width W' for the 10 standard clusters against [Fe/H] values 
from lKraft Ivansl ( 2003 ) derived using MARCS model atmospheres. The fitted line has a slope 
of 0.549 dex/A. 
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the values of W8542, Ws662 and their sum, together with their associated uncertainties. 



3.2. M22 



For the 55 stars in the M22 field with reduced spectra, 51 have velocities compatible with 
cluster membership. The mean velocity of these stars i s -150 ± 2 km s" 1 which agrees well with 
the value -148.9 ± 0.4 km s _1 tabulated by lHarria (199a). The velocities of the other 4 stars clearly 
exclude them as members. The W§542 + ^8662 values for the member stars were then determined 
in the same way as for the standard cluster stars. A plot of the values against V — Vhb is shown 
in Fig. [H Shown also on the Figure is the best fit of a line with a slope of -0.51 A/mag; as for 
the standard clusters only the 41 M22 members with V — Vhb < 0.2 were included in the fit. For 
completeness we note that a least squares fit to these 41 points yields a slope of -0.49 ± 0.06 A/mag 
fully consistent with the slope derived from the standard cluster observations. 

The average W value for the 41 M22 stars with V-V H b < 0.2 is 2.912 ± 0.044 (std error of the 
mean), which corresponds to an abundance \Fe/H]MAacs = -1-77 dex using the calibration given 
in equation (1) above. The formal uncertainty in this value is small (0.549 x 0.044 = 0.02 dex), but 
in practice the actual uncertainty will be of order of the calibration uncertainty, i.e. approximately 



0.1 dex . Comparison of this value with other determinations is not straightfo rward. iRutledge et al 



( 199 7_a|) did not obse r ve M2 2 and so there is no value for this cluster given in iKraft &: ivansl (12003). 



Carretta & Grattonl (| 19971 1 give [Fe/H] = -1.48 ± 0.06 on their scale but iKraft & Ivansl (|2003h 
note that using MAR CS models, they der i ve ab undances that are syste matically low e r by a bout 
0.2 dex than those of ICarretta fc Grattonl (119971 ) for the same clusters. IZinn West! (119841 ) give 
[Fe/ H] = -1.75 ± . 08, w hich is in close accord with our determination. This is not surprising 



since 



Kraft &: Ivansl (120031 ) note that at metallicities of this order, the diff erence betwee n thei r 



MARCS model abundances (whose scale we have adopted here) and those of IZinn k, Westl (119841 ) 



are small. Ilvans et al 
scale, while 



(1 2 0041 ) uses an abundance of [Fe/H] = -1.7 dex, presumably on the MARCS 
Harrisl (|l996^ = —1.64 for M22. These are again consistent with our 

determination. Similarly iMarino et al.l (|2009l ) give a mean metallicity for M22 of [Fe/H] = -1.76 ± 



0.02 (internal errors only, weighted mean of the UVES and GIRAFFE data), essentially identical 
to our determination. 

We give in Table [2] the ID numbers and positions for the 51 M22 member stars, along with 
the measured heliocentric radial velocities in km s _1 , the magnitude difference from the horizontal 
branch, AV, the individual Ws542 an d Wse62 measurements in A and their errors, as well as their 
sum, X W, and its uncertainty. 
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Table 1. Standard Cluster Data 



ID 


RA (2000) 


Dec (2000) 


Rad Vel 


AV 


W 85 42 


e 


Wg662 


e 


EW 


e 








(km s -1 ) 


(mag) 


(A) 




(A) 




(A) 




M4 1_286 


16 23 43.2 


-26 31 54 


87 


-0.53 


2.57 


0.15 


1.87 


0.09 


4.44 


0.17 


M4 1.1227 


16 23 41.0 


-26 31 30 


69 


-0.58 


2.55 


0.12 


1.94 


0.09 


4.49 


0.15 


M4 1.1604 


16 23 54.8 


-26 31 20 


53 


-0.23 


2.31 


0.12 


1.81 


0.09 


4.12 


0.15 


M4 1.3266 


16 23 38.6 


-26 30 38 


58 


-0.41 


2.52 


0.17 


1.77 


0.11 


4.29 


0.20 


M4 1.4341 


16 23 41.0 


-26 30 09 


81 


0.21 


2.41 


0.19 


1.57 


0.11 


3.98 


0.22 



Note. — This table is available in its entirety in a machine-readable form in the online journal. A 
portion is shown here for guidance regarding its form and content. 



Table 2. M22 Member Data 



ID 


RA (2000) 


Dec (2000) 


Rad Vel 


AV 


W8542 


e 


W8662 


e 


sw 


e 








(km s _1 ) 


(mag) 


(A) 




(A) 




(A) 




1.232 


18 36 41.5 


-23 54 39 


-148 


-0.10 


1.68 


0.06 


1.47 


0.06 


3.15 


0.08 


1.756 


18 36 41.1 


-23 54 33 


-136 


0.98 


1.46 


0.15 


0.97 


0.11 


2.43 


0.19 


1.1211 


18 36 42.3 


-23 54 28 


-144 


0.52 


1.73 


0.11 


1.30 


0.07 


3.03 


0.13 


1.2276 


18 36 37.4 


-23 54 17 


-146 


-0.42 


1.61 


0.07 


1.32 


0.05 


2.93 


0.09 


1.3666 


18 36 36.3 


-23 54 02 


-148 


-1.57 


2.00 


0.09 


1.55 


0.05 


3.55 


0.10 



Note. — This table is available in its entirety in a machine-readable form in the online journal. 
A portion is shown here for guidance regarding its form and content. 
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M22 




1 0-1-2 -3 

v - v HB 

Fig. 4. — Plot of Ca II line strength (M / 8542+W / 8662) against magnitude difference from the horizontal 
branch (V — Vhb) for the 51 M22 members. The solid lines, in order of increasing Ws542 + Ws662 
values at V — Vhb = -1.5, are the relations for the standard clusters M15 ([Fe/H] marcs = ~2-45), 
NGC 4372 (-2.29), NGC 6397 (-2.12), M10 (-1.48), M5 (-1.32) and M71 (-0.82). The dot-dash 
line is the fit of a line of slope -0.5lA/mag to the M22 stars with V — Vhb < 0.2. Vertical bars on 
each point show the measurement uncertainty in the line strength. 
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4. Discussion 

The most striking thing about the distribution of points in Fig. 2] is the large apparent scatter 
about the fitted line: the rms dispersion in W§542 + ^8662 at fixed V — Vhb for the 41 stars with 
V — Vhb < 0.2 is 0.28 A, which is considerably larger than the mean measurement error, 0.11 A. 
This stands in stark contrast to the situation for the standard clusters. For example, the stars 
observed in the standard clusters M15, NGC 6397 and M10 have mean measurement errors of 0.11, 
0.11 and 0.14 A, respectively, values comparable to that for M22. Yet for these clusters the rms 
dispersions about the fitted line are 0.16, 0.13, and 0.08 A, comparable to errors and substantially 
less than is the case for M22. It is also clear from Fig. E]that the fainter stars in M22 show a similar 
dispersion. We now investigate possible origins for this large spread. 

The first possibility is that the sample is contaminated with non-members. This seems very 
unlikely: the observed velocity dispersion of the stars classified as members is 12 km s" 1 , which 
is similar to the observed dispersion for other member sets. For example, the velocity dispersions 
of the member stars in the standard clusters M15, NGC 6397 and M10 are 20, 12 and 8 km s _1 , 
respectively. Further, the M22 stars whose measured velocities lie furthest from the mean are not 
distinguished from the remainder of the sample in Fig. SJ nor do the stars that stand furthest from 
the mean line in the figure have velocities that are notably discrepant. We conclude that all the 
stars in Fig. [J] are indeed likely members of M22. 

The second possibility is that the large scatter is produced by differential reddening across 
the observed field. Such reddening variations would induce scatter in the V — Vhb values and 
thus potentially increase the scatter about the mean line. To investigate this possibility we first 
examined the spatial location of the observed stars. Specifically, we have split the sample into two 
groups, those that lie above the mean line (potentially more reddened on average, 19 stars) and 
those that lie below the mean line (potentially less reddened, 22 stars). Figure [5] shows the outcome: 
the first group are plotted as red stars while the second group are plotted as blue squares. There 
is clearly no straightforward segregation of the red and blue points in the Figure. Consequently, if 
reddening variations are the explanation, then the variations must occur on small scales of order 
~20-30", perhaps less. 

While it is unclear whethe r reddening variations on such scales are present in our M22 field, 



we note that lLvons at al.l (|1995l ) have shown that E(S — V) variations of order 0.05 mag can occur 
on scales as small as few arcsecs in the field of the globular cluster M4. Thus we need to at least 
enquire into the size of the reddening variations required to produce the scatter in Fig. |H We have 
done this via Monte-Carlo simulations. We start with the assumption that there is no intrinsic line 
strength spread other than that induced by the measurement errors in the Ws542 + W / 8662 values. 
The trials are conducted as follows: for each of the 41 V — Vhb values with V — Vhb < 0.2, 
the mean line in Fig. |4] is used to generate initial Ws542 + ^8662 values. These initial values are 
then altered by random draws from a gaussian distribution with mean zero and standard deviation 
equal to mean error in the observed Wg542 + Ws662 values, which is 0.11 A. In the absence of any 
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additional broadening mechanism, the average over many trials would have a dispersion about a 
fitted line of slope -0.51 A/pix equal to the mean observational errors, and a W value equal to 
the observed value. We verified that this is indeed the case. To simulate the effects of differential 
reddening we then perturbed the V — Vhb values with random draws again from a gaussian, keeping 
the Wg542 + W / 8662 values (after perturbation with the error distribution) unaltered. This second 
gaussian has mean zero and a standard deviation as input, representing a(E(B — V)). After the 
V — Vhb values are perturbed, assuming Ay = 3.2 E(B — V), a line of slope -0.51 A/pix is fit 
taking care to exclude any stars whose perturbed V — Vhb value exceeds 0.2, and the corresponding 
W' and dispersion about the line calculated. We find that in order to consistently reproduce the 
observed dispersion of 0.28 A with mean measurement errors of 0.11 A, we require a(E(B — V)) ~ 
0.12-0.15 mag. Clearly adopting a larger value than 3.2 for R = Ay /~E(B — V) would increase the 
dispersion about the mea n line for a given a (~E(B — V)). However, we find that even adopting R 
= 3.6 (cf. Appendix F of lBessell et al.lll998l ). the result is not materially altered — the minimum 
a(E(B — V)) required for the observed dispersion to be consistently realized in the trials becomes 
0.11 rather than 0.12 mag. 

These values of a(E(B — V)) are considerably lar ger than existing esti mates of the degree of 
differential reddening in the field of M22. For example, iMonaco et all (J200J) estimate a maximum 
range in K(B — V) across the entire face of the cluster as 0.06 mag, larger variations would be 
inconsistent with their data. We conclude therefore that there is an intrinsic spread in the Ws542 + 
^8662 values for M22, over and above that due to any differential reddening. 



4.1. The Abundance Spread in M22 



If the intrinsic spread in the 1^8542 + ^8662 values for the M22 stars is not due to non- 
member contamination and is not primarily the result of differential reddening, then the remaining 
alternative is that it is the consequence of an intrinsic abundance spread in the cluster. Noting that 
for each individual star we can define a W value as W8542 + H"8662 + 0.51(V — Vhb), equation (1) 
can then provide an abundance estimate [Fe/H] MARCS for each individual star. We show in Fig. 
[6] a generalized histogram made from the individual [Fe/H] a/ arcs values and the corresponding 
abundance uncertainties that follow from the measurement uncertainties in the Ws542 + Ws662 
values. These have a mean value of 0.06 dex. We also show in the figure the contributions to the 
total from the stars below (more metal-poor) and above (more metal-rich) the fitted line in Fig. 
IH which corresponds to an abundance [Fe/H] = -1.77 dex. The method assumes that there is no 
significant contribution to the abundance determinations from differential reddening. This is likely 
to be the case fo r the following reason. If the differential reddening has a range AE(B — V) ~ 0.06 



- 0.08 mag (e.g.. IMonaco et alj|2004l ; lAnthonv-Twarog et al.lll995l ). then presumably a(E(B — V)) 
~ 0.02 and thus a (Ay) ~ 0.06 mag with the same value for a(V — Vhb)- This converts to <r([Fe/H]) 
~ 0.02 dex, which is much smaller than the observed value of <r([Fe/H]) = 0.15 dex. This contrasts 
with the photometric case where a(V — I) from differential reddening of size a(E(B — V)) ~ 0.02 
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Fig. 5. — Plot of the positions of the 19 M22 member stars that lie above the mean relation in 
Fig. S] (red stars) and of the 22 stars that lie below the relation (blue squares). The 6.8' x 6.8' 
field-of-view from which the observed stars were selected is outlined by the dotted lines. The cluster 
center is shown by the green circle. Note that there is no obvious spatial separation of the red and 
blue symbols. 
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mag, and a(V — I) from cr([Fe/H])j n t ~ 0.15 dex are comparable in size and thus more difficult to 
distinguish. 

Figure [6] clearly shows that there are multiple components to the abundance distribution in 
M22 — it rises rapidly to a narrow peak at [Fe/H] a/ arcs = -1-88 yet there is a broader tail towards 
higher abundances. We can quantify this distribution in a number of ways. For example, for the 
22 stars in metal-poor group, the inter-quartile range is only 0.05 dex, while for the 19 stars in the 
metal-rich group, the inter-quartile range is notably larger, 0.15 dex. The median abundance for 
this group is [Fe/H] marcs = -1-64 and the most metal-rich star has [Fe/H] m ARCS = -1-43 ± 0.07. 
The most metal-poor star has [Fe/H] marcs = — 2.1 ± 0.05 and it is 0.14 ± 0.07 dex more metal- 
poor than the next most metal-poor star. Whether this star represents a separate third metallicity 
grouping at low abundances and with a small fraction of the total (few percent at most), or is 
simply a statistical outlier, cannot be determined without a larger sample of member stars. For the 
entire sample the inter-quartile range is 0.24 d ex, which is comparable to the M22 abundance range s 



found by the earlier spectroscopic studies of iNorris fe Freeman! (119831 ) and iLehnert et al 



though it is clearly in conflict with the results of llvans et al.l (12004! ). Ou r results are summa rised 



(1991 



in Table [3l We defer to the next section a comparison with the results of lMarino et all (|2009l ) . 



We can also take 'toy' models for the intrinsic abundance distribution and compare them 
with the observed distribution after convolution with the measurement errors. One such model 
is shown in the upper panel of Fig. [7l Here we have assumed that the M22 stars consist of two 
populations: one with abundances uniformly distributed between [Fe/H] = -1.95 and -1.83 making 
up 44% of the total, and a second with abundances uniformly distributed between -1.83 and -1.50 
dex and 56% of the total. This abundance distribution is shown in the insert in the upper panel 
of the figure. Convolution with the measurement errors then gives the dot-dash line, which is 
an acceptable representation of the observed data. The contributions of the two components are 
shown by the dotted lines. We have not conducted an exhaustive parameter search so it is likely 
that other intrinsic abundance distributions could produce similar fits to the observations, although 
two components with different abundance ranges does seem to be a requirement. We note that the 



Table 3. M22 Abundance Data 



Sample 


N 


([Fe/H]) 


a obs ([Fe/R}) 


IQR 


All 


41 


-1.77 


0.15 


0.24 


Metal-Rich Group a 


19 


-1.63 


0.09 


0.15 


Metal-Poor Group b 


22 


-1.89 


0.07 


0.05 


a Stars with [Fe/H] 


> - 


1.77 






b Stars with [Fe/H] 


< - 


-1.77 







-2.2 -2 -1.8 -1.6 -1.4 -1.2 

[ Fe / H Ws 

Fig. 6. — A generalized histogram (solid line) of the individual [Fe/H] m a rcs values for the 41 M22 
stars with V — Vhb < 0.2. Each value has the abundance uncertainty from the uncertainty in the 
W8542 + W-8662 measurement. Shown also are the generalized histograms for the 22 stars below the 
fitted line in Fig. |4] (dotted line), and for the 19 stars above the fitted line (dot-dash line). 
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model does not reproduce the extreme metal-poor tail of the observed distribution adequately. A 
third component could no doubt be added to the model to fix this, but given that the metal-poor 
extension results from a single star, the addition of a third component does not seem warranted at 
this stage. As noted above, a larger sample of member stars is needed to fully characterize the most 
metal-poor part of the distribution. We note also that assuming a single abundance for the metal- 
poor group does not fit the observations adequately — a range in abundance in this population 
is apparently required. However, we have not included the effects of differential reddening, which, 
as outlined above, could induce cr([Fe/H]) ~ 0.02 dex. Consequently, the extent of the metallicity 
range in the metal-poor group might well be smaller than assumed in this toy model. 

We can also compare the observed abundance distribution with the predictions of simple 
models of chemical evolution. One of the simplest such model is that in which star formation 
proceeds in gas whose initial abundance is Zq under the assumption of instantaneous recycling 
and with the rate of gas loss from the system proportional to the star formation rate. In such 
a model the metallic ity distribution function f(Z) is ch aracterized solely by Zq and by the mean 
abundance (Z) (e.g., lHartwicid 119761 : iNorris et al.lll996l ). In the lower panel of Fig. [7] we show a 
comparison of the prediction of such a model with the M22 observations. The model has been 
calculated with log(Zo/Z sun ) = -1.95 and log((Z)/Z sun ) = -1.77 dex, which is the observed sample 
mean abundance. We assume Z sun — 0.017. The model distribution was then convolved with the 
0.06 dex mean [Fe/H] error. The fit to the observations is adequate, but it does seem likely that 
to improve the fit , a second component with higher values for Zq and (Z) would be required (cf. 
Norris et al.lll996l ). However, it is not clear what additional insight would result from computing 
such a two component model given that it does not help understand how such a model might arise 
physically. 

Additional insight can be gai ned by comparing the M22 abundance distribution with that 



Norris et al 



for u Cen (e.g., INorris et al.lll996 ) . To do this we have taken the [Ca/H] abundances for ~500 



(1996) and generated a generalized histogram using the average 



lo Cen red giants from 

abundance error given by INorris et al.l (|1996l ). which, at 0.05 dex, is very similar to that for the 
M22 observations. The u Cen histogram was th en shifted first by -0.4 dex , which corresponds to 
the mean [Ca/Fe] ratio for uj Cen red giants (e.g. . INorris &: Da Costalll995bl ). and then by a further 
-0.09 dex to match the [Fe/H] of the peak of the M22 abundance distribution. It was then scaled 
so that the peak height coincided with that for M22. The outcome is shown in Fig. El The two 
distributions are not significantly different on the metal-poor side of the peak given the smaller M22 
sample. It is clear, however, that the u> Cen distribution is broader on the metal-rich side, though 
the 'rate of decline' away from the peak is similar. It also appears that the 'second abundance 
peak' in M22 is closer in abundance to the main peak than is the case for uj Cen. Furthermore, 
the u Cen abundance distribution continues to significantly higher abundances than M22, a result 
that is not obviously an outcome of the smaller size of the M22 observed sample. Despite these 
differences, the general similarity between the two distributions suggests that they could readily be 
the result of the same physical process or processes, except that it has gone on longer in uj Cen, or 
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Fig. 7. — (a,) Upper panel. The solid line is the generalized histogram of the individual M22 
[Fe/H]MARCS abundances from Fig. [6j The dot-dash line, which is a reasonable representation of 
the observations, is the result of convolving the abundance distribution shown in the insert with 
the measurement errors. This distribution has two components, one with abundances distributed 
uniformly between [Fe/H] marcs = -1-95 and -1.83, and the other distributed uniformly between 
-1.83 and -1.50 dex. The population ratio is 1 to 1.25. The contributions of each component are 
shown as dotted lines, (b) Lower panel. The solid line is again the generalized histogram of the 
individual M22 [Fe/H] marcs abundances from Fig.[6l The dot-dash line is the abundance distribu- 
tion for a simple chemical model in which the rate of gas loss is proportional to the star formation 
rate. The model assumes an initial abundance of log(Zo/Z SMn ) = -1.95 and has log((Z)/Z sun ) = 
-1.77 dex corresponding to the mean abundance of the observed sample. 
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terminated sooner in M22, perhaps because of the difference in mass between the systems. 



Given that, in contrast to M22, most globular clusters are homogeneous with respect to heavy 
element abundances, the similarity between the M22 and uj Cen abundance distributions seen in 
Fig. [8] naturally leads to the following speculation. The unusual properties of u Cen have frequently 
led to the i nterpretation th at the stellar system is the nuclear remnant of a now disrupted dwarf 
galaxy (e.g., Freeman 19931 ). Consequently, it is plausible to suggest that M22 may also be a nuclear 
remnant of a disrupted dwarf galaxy. The fact that M54, the central star cluster of the Sgr dwarf, 
also shows an internal abundance range lends further support to the suggestion. The stellar system 
lu Ce n lies in a tightly bo und r etrograde orbit wh i ch nev er rises very far above the Galactic Plane 
(e.g., iDinescu et al.lll999l ). Yet iBekki Freeman! (|2003l ) have shown that no inconsistencies arise 
in evolving a nucleated dwarf galaxy in the tidal gravitational field of the Galaxy from an initial 
large galacto centric distance to a nuclear r emnant with an orbi t resembling that of the present-day 
lj Cen. As regards M22, t he compilation of Dinescu et al. (jl999l ) shows that its orbit is fairly typical 
for inner halo objects (cf. lCarollo et al.1120071 ). It is strongly prograde (G = 178 ± 20 km s _1 ) with 
apo- and pericentric distances of a pproximately 9.5 an d 2.9 kpc, respectively, and a maximum 
height above the plane of 1.8 kpc (jDinescu et al.lll999l ). Given the results of IBekki &: Freeman 
([2003) for uj Cen, it seems likely that one could also start with a nucleated dwarf galaxy at large 
galactocentric distances and evolve it such that the nuclear remnant has an orbit similar to that of 
M22. 



It is also noteworthy that both w Cen and M22 have strong blue horizontal branch populations, 
as does M54. In the case of M22 this is not unusual given its relatively low mean metal abund ance 
and relatively s mall distance from the Galactic Center: in the terminology of IZinnl (| 19931 ) and 
Lee et al.l (| 19941 ) M22 is an 'old-halo' cluster. But it is also interesting to note that M22, u Cen and 
M54, the three clusters for which there is definite ev i dence f or the presence o f internal abundance 
ranges, are all in the 'EHB' classification of lLee et al.l (|2007l ). iLee et al.l (120071 ) have shown that the 
EHB clusters form a kinematically distinct group, and argue that such clusters have their origin 
as the cores or central star clusters of stellar systems accreted during the earliest phases of the 
formation of the Galaxy. This is then consistent with our interpretation of M22 and its unusual 
properties. 



4.2. Comparison with the results of iMarino et al.l ([2009) 



In a recent paper IMarino et al.l (|2009l ) have analyzed high dispersion UVES spectra for 17 red 
giants in M22, supplemented by lower resolution GIRAFFE spectra for a further 14 stars (one 
of which is in common with the UVES set). The results can be summarized as follows. The 
distribution of the abundances of the s-process elements Y, Zr and Ba in the UVES sample is 
apparently bimodal, with one group of stars showing significantly larger average values of [(Y, Zr 
and Ba)/Fe] compared to the other. Specifically, in the UVES sample, 7 of the 17 stars (~40%) 
have ([s/Fe]) = 0.46 ± 0.03 while the remainder have ([s/Fe]) = 0.02 ± 0.02 dex (errors are internal 
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[Fe/H], 



Fig. 8. — The solid line is the generalized histogram of the individual M22 [Fe/H]^/ arcs abundances 
from Fig. [6l The d ot-dash line is the abundance distribution generalised histogram for to Cen from 
Nor r is et al.1 ( 1996 ). It has been shifted horizontally and scaled vertically to the peak of the M22 
distribution. 



-24- 



only). The [Ba/Fe] values for the stars observed with GIRAFFE are consistent with this result. 
Moreover, for the UVES sample, the s-process rich group has both higher average [Fe/H] and [Ca/H] 
abundances, by 0.14 ± 0.03 dex and 0.25 ± 0.04 dex, respectively. For the GIRAFFE sample the 
difference in the average [Fe/H] values is 0.18 ± 0.04 dex. The interquartile range in [Fe/H] for the 
full sample is 0.16 dex, with the most metal-poor star at [Fe/H] = -1.94 and the most metal-rich 
at [Fe/H] = - 1.59, with (interna l ) unce rtai nties of 0.09 dex. Th ese values are reminiscent of the 
earlier work of lNorris &; Freeman! (|1983l ) and lLehnert et al.1 (|199ll ). and are entirely consistent with 
the results presented here. 



There are 3 stars in common between the 



Marino et al.l (|2009l ) UVES observations and the 



sample presented here, with iMarino et al.l (|2009l ) stars 200068, 200101 and 200083 corresponding 
to our stars 1 16051, 1_3683 and 1 ,10101. For these stars, which have [Fe/H] values of -1.84, -1.7 4 



and -1.63 in 



Marino et al.l (120091 ) , the differences in [Fe /H] values in the sense of IMarino et al 
|200j|) abundance minus that derived here, are 0.09 ± 0.11, 0.02 ± 0.13 and 0.01 ± 0.11 dex, 



where the error given is the combined internal uncertainty. Clearly, although the sample is small, 
there is no evidence for any systematic offset (the mean difference is 0.04 dex) and the standard 
deviation of the differences, ±0.05 dex, is consistent with the errors. Consequently, we have used 



a (two -sided) K-S test to compare the [Fe/H] distribution for the 30 stars in the IMarino et al 



(|2009l ) full sample with that for the 41 stars with [Fe/H] values measured here. We find that 
the null hypothesis, namely that the two samples come from the same underlying distribution, 
cannot be ruled out, even at a 10% sign ificance level. This has one immediate consequence. As 
noted in the Introduction, iPiottd (J2009J) has shown that in a ffi>T-based c-m diagram the M22 
sub-giant branch has a bimodal structur e with the upper (b righter) branch contributing 38 ± 5% 
and the lower 62 ± 5% of the total (see IMarino et al.l |2009j) . One interpretation of this bimodal 
structure in the c-m diagram is that it reflects the abundance distribution in the cluster, with the 
fainter sub -giant branch corres ponding to the more metal-rich population, with a negligible age 
difference. IMarino et al.l (|2009l ) note, however, that for a difference in average [Fe/H] of 0.14 dex, 
the predicted difference in the sub-giant branch lo cation is ~0.1 mag in F606W, smaller than 
the observed difference of ~0.17 mag. Further, in IMarino et al.l (|2009l ). the metal-rich group is 
less than half the sample (37 ± 13%), in contrast to the situation in the c-m diagram. Our data 
alleviate both of these concerns. For example, the difference in average abundance between the 
two populations shown in the insert in the upper panel of Fig. [7] is 0.23 dex, which would predict 
a bigger magnitude difference on the sub-giant branch. Also, the ratio of the populations, 56% 
metal-rich to 44% metal-poor, is more in accord with the observed ratios in the c-m diagram. 
Detailed modeling of the c-m diagram is required to place limits on the abundance distribution, 
and variations in other plausible parameters such as total CNO and age, required to reproduce the 
observations. 

One further point can be made concerning M22 and u Cen. We argued above that the over- 
all M22 abundance distribution is similar to that of uj Cen, except that M22 lacks the tail to 
higher abundances. The availability of high dispersion abundance analyses for M22 stars from 
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Marino et al.l (j2009l ) means that the comparison can n ow be extended to i ndivid ual abundance 



ratios, using u Cen red giant results from, for example, iNorris Da Costal (|1995bl ). We show an 
exam ple in Fig. [9l The upper pa nel shows [Ba/Fe] as a function of [Fe/H] for 40 oj Cen red giants 
from INorris fe Da Costal (|1995bl ). The [Ba/Fe] ratio rises with increasing [Fe/H] until [Fe/H] ~ 
—1.3 after which the Ba and Fe abundances change in lock-step. The dot-dash line is a least squares 
fit to the stars with [Fe/H] between -1.9 and —1.3, e xcluding the three metal-p oor stars with high 
[Ba/Fe], at least one of which is a carbon-star (see INorris feDa Costal Il995bl ). The lower panel 
shows the [Ba/Fe] ratios for the 30 M22 stars in iMarino et al.1 ^OOfll ). The dot-dash line is repro- 
duced from the upper panel except that it has been shifted to lower [Fe/H] by 0.2 dex. While there 
are differences in detail, such as the metal-poor stars in M22 having [Ba/Fe] ~ 0.0 dex but -0.2 
in Cen, the relations are clearly similar, indicating that it is likely similar enrichment processes 
went on in both clusters. 



5. Conclusions 

We have shown here that there is strong evidence for the existence of a significant i nterna l 



abundance range in the globular cluster M22, consistent with the results of IMarino et al.l ([2009). 
The cluster thus joins u Cen and M54 as the only Galactic globular clusters in which intrinsic [Fe/H] 
ranges have been established. The M22 abundance distribution rises sharply from low abundances 
to a distinct peak at [Fe/] ~ -1.9, with a broad tail to higher abundances. Consequently, it is 
probable that at least two components are needed to describe the distribution. The abundance 
distribution also bears a qualitative similarity to that for to Cen, although the overall scale in M22 
is considerably smaller. Nevertheless, the similarity suggests a comparable origin for the abundance 
ranges in both systems. Indeed, given that M54 is currently the central star cluster of the Sgr dSph, 
and that oj Cen is frequently postulated as being the remnant nucleus of a disrupted dwarf galaxy, 
it seems logical to conclude that M22 may also be the remnant nucleus or nuclear star cluster of a 
disrupted dwarf galaxy. 
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